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A generic saddle-node bifurcation is proposed to modelize fast transitions of finite amplitude arising in geo- 
physical (and perhaps other) contexts, when they result from the intrinsic dynamics of the system. The fast 
transition is generically preceded by a precursor phase which is less rapid, that we characterize. In this model, if 
an external source of noise exist, the correlation length of the fluctuations increases before the transition, and its 
spectrum tends to drift towards lower frequencies. This change in the fluctuations could be a way of detecting 
catastrophic events before they happen. 

PACS numbers: 



I. STATEMENT OF THE PROBLEM 

Earthquakes, like volcanic eraptions as well as other physical phenomena and perhaps also some kinds of socio-economical 
"revolutions", show an abrupt transition from one state to another. We consider cases where this transition is intrinsic (not the 
result of an excitation from outside) and dynamical in the sense that, as a parameter changes slowly, the system jumps by a 
finite amount in a time much shorter than the typical time of evolution of the external parameter In earthquake physics, this 
typical time of evolution, the earthquake recurrence time, is on geological scales of plate tectonics although the time scale of 
seismic ruptures is within the second to minute range [1]. Our basic assumption is that, as a dynamical system, an earthquake 
shows a "dynamical saddle-node" bifurcation. At the bifurcation point, a pair of fixed points, one locally stable the other locally 
unstable, merge and vanish as a control parameter varies. Take a damped dynamical system, with a coordinate x{t) solution of 
the equation 

^ = (1) 
dt dx 

In this equation V{x) is a potential. In the geophysical context of earthquakes, the scalar variable x could be the relative 
displacement across the fault. The equation ([T]i is too general to be very helpful. However, as time varies slowly, it may describe 
a saddle-node bifurcation where a stable equilibrium disappears, assuming that V depends slowly on time in a prescribed way, 
to become a function V{x, t). Near the transition, one may use a mathematical picture which is correct for a short time around 
the transition if the potential V{x, t) is a smooth function (see below for what happens beyond this local study). 
Assume first that V{.) does not depend explicitely on time and takes the form 

V{x)=-{^x^ + bx), (2) 

with h real constant (for the moment). 

For b negative V{x) has two real extrema {i.e. the roots of ^ = 0), one —^/~-b is a stable equilibrium, the other, ^/~-b, is an 



dx 

unstable equilibrium. For 6 = the two equilibria merge and disappear for b positive, see Figure ([T}a). This is the saddle-node 
bifurcation. The shape of V{x) near a; = and for b small is universal: for a given smooth V{x) showing this saddle-node 
bifurcation, one can always rescale time and external parameter to find the "local' problem in this form. 

The extension to time dependent variation of the control parameter b goes as follows. If 5 is a smooth function of time, one 
can assume that b{t) crosses the critical value, i.e. zero in the present case, at time zero in such a way that b{t) = at + ... with a 
non zero constant, and the dots being for higher terms in the Taylor expansion of b{t). For t and x close to zero, after elementary 
rescaling, one can represent the dynamical system ([1}, close to the saddle-node bifurcation, by an "universal" parameterless 
equation 



-=x^+t. (3) 

Outside of the neighborhood of a; = 0, the solution of ([T]i depends on other parameters like the one defining V{.) far from 
X = , as studied below. 
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FIG. 1: (a) Cubic potential for b = -1, 0, 1. (b) Quartic potential, b = -1, 0, 1 



A. Solution close to the saddle-node 



Let first consider what happens close to the saddle-node bifurcation. We look for a solution of equation (|3) transiting from the 
"stable" fixed point at "large" negative times to the rolling down one towards positive value of x at positive times. This solution 
behaves like x{t) w —\f~-i at large negative times. The equation (O is of the Riccatti type and can be integrated by introducing 
the function y(i) such that x(t) = — — where y' = ^ and y{t) is a solution of Airy's equation. 



y" + ty = 0. (4) 
The solution relevant with the given condition for x{t) att ^ — oo, is the Airy function Ai{—t) which writes 

Y{t) = Ai{-t) = / cos(— - ut)du, (5) 
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and leads to the curve x{t) drawn on Figure(|2]-a). Yet we have only solved the transient problem near the saddle-node bifurcation. 
The transition ends-up when t becomes equal to the first zero of the Airy function Ai{—t), corresponding to a divergence of the 
original x{t). Let tc be this critical value of t, i.e. the smallest root of Y{t) = 0, a pure number, about tc ~ 2.338, and let us 
look at the behaviour of x{t) just before this transition. From the Laurent expansion of Y{t) close to tc, and returning to the 
original variable x{t), one obtains 

a^W«^-|(ic-0-- (6) 

As this solution diverges, it looses its validity because the "universal" dynamical equation (|3]l was derived under the condition 
that X remains close to zero. This local theory cannot deal with finite variations away from the critical conditions, therefore we 
shall add finite amplitude effects to limit the growth of the instability after the transition. 

We shall study now two questions, first the response of this dynamical system to an external noise source, then the dynamics 
of a system showing a saddle-node bifurcation of the type just studied and reaching a new stable fixed point after this bifurcation. 



B. Response to an external noise source 

We explore first the response of our system to a small external noise, and search whether the response to the noise changes 
qualitatively and so could be a signal ahead of the transition. 

Let us consider the equation Q with a small noise added, so that equation (|3]l is replaced by 



^=x2+i + eC(i), (7) 

where C,{t) is a random function of time, and e a small factor 

In the limit e small, one can solve equation (|7]i by expansion in powers of e, 

x{t) = xa{t) + exi{t) + .... (8) 

where 

Y'(t') 

Mt) - -Y(tf, (9) 
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The linear response to the noise is 



xiit) 



(10) 



Because Y'^{t) tends rapidly to zero as (t) tends to minus infinity, one can take to = — oo to get rid of the effect of the initial 
conditions. 

Let us take a delta-correlated (or white) noise, such that 



(C(ia)C(4)> = S{ta - h). 

Note that, including with e = 1, the solution of equation O with noise is still very close to the noiseless solution xo(t). 
The correlation function of xi{t) is given by 



(11) 



{xi{t)xi{t')) 



al{t,t') 



(12) 



whose behavior for large negative values of both t and t', is derived from the asymptotic expression of the Airy function 
Ai{—t) « 2^y^{-t)^/'i • 'Ss'^'^iiig ^ ~ J' ^'^'1 F{'w) = 1 — w"^^^, the variance = {{x{t) — xo(t))^) of the fluctuations writes 



{x^(tf)^{-t) 



In the limit {—t) oo the integral is concentrated near w = 1 so that 



(13) 



(14) 
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(b) 




FIG. 2: (a) Solutions of equation ([7]l, with and without noise, e = (smooth curve) and e — 1 (noisy curve), (b) Noise (^(t). (c) Solution of 
equation l il7t for a — 10"^. The rectangle around the origin defines the region —to < t < to and — 1/to < x < 1/io, with to ~ a~^^''^. The 
critical time is tc ~ 2.34 to- The two vertical lines inserted between the two aiTows delimitate the large slope time duration, of order unity. 



which shows that the fluctuations increase as time goes on, some time before the transition itself. As the transition approaches, 
the variance of the fluctuations increases close to the critical time tc, because Y{tc) ~ 0. 

Because of the divergence of the solution at t = tc, it does not make sense to describe the dynamical behavior of the 
fluctuations due to the external noise very close to tc- This unbounded growth is a consequence of the local cubic form of V{x) 
as expanded near x = 0, which is valid around a; = only, in obvious contradiction with the fact that x{t) diverges. 



II. DYNAMICAL SADDLE-NODE BIFURCATION TOWARD A NEW STEADY STATE 



To suppress the divergence of x{t) after the saddle-node bifurcation we add a stabilizing (positive) term to the potential V{x) 
which becomes quartic. 



x^ 



V,ix) = -—-bx+—, (15) 
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as drawn in Figure([T}b). Because of the growth of Vq{x) at infinity, like x^, the solution of the differential equation 

does not diverge at finite time. The equation (flSl l can be written in the given scaled form, provided the coefficient of is 
positive. For such a potential one parameter only remains. In equation ( fTSl ) the coefficient is chosen as b, the one of the linear 
term. For 6 = the dynamical system ( fTSl ) is exactly at the saddle-node bifurcation, because at b = x = the first and second 
derivative of Vq{x) both vanish, but not the third one. Contrary to the case of the pure cubic potential, this system has always, 
that is for any value of b, a stable fixed point beyond the pair of fixed points collapsing at the saddle-node bifurcation. This 
makes it a fair candidate for describing the dynamical saddle-node bifurcation without blow-up. 

As in the previous case, we shall look now at the case of a time dependent b, that will be taken as b = at with a positive 
constant. Because of the rescaling of the cubic and quartic term, the parameter a cannot be eliminated (another possibility is 

3 4 

to put a parameter in front of the cubic term). For the potential Vq{x) = — atx + ^ we shall analyse the solution of the 
dynamical equation 

^ = at + x'^-x^, (17) 
at 

tending at large times to the quasi-equilibrium point x = (at)^/"^, t being considered as a parameter, see Figure(|2}c). Moreover 
we consider the limit a small, which could describe a wide range of slip phenomena ||2tl, as earthquakes where a is generally 
very small, of order 10~^ (see below). 



A. Three time ranges for small a 

In this limit we show first that there are three characteristic time intervals, depending how x is close to zero. 

The long time scale is the average recurrence time of an earthquake at the same site along the fault. It is typically of order 
^phys ^ years. In our model it is the time needed for the potential Vq{x, t) to change significantly, to move from a pair of 
fixed points to a saddle-node bifurcation. Because time enters in Vq{x, t) through the combination {at), the adimensional time 
needed for a change of shape of Vq is of order 

tb^a'\ (18) 

The short time icqk is of order unity in our model equation (fTTt as stated in the next paragraph. It is the duration of the abrupt 
change of the slope of the solution x{t). This short time corresponds to the dynamic rupture duration, which is typically of order 

ten seconds for a magnitude 6 earthquake, t^^ ~ IDs 111]. Therefore the ratio of these two time scales = -j^, is small as 

a — 1.6 10^^ in the geophysical context. 

There is another time scale, to, the time interval standing before the transition, and close to it, during which the potential is 
very flat, while the solution has not jumped. During this time, x and at are much smaller than unity, then the cubic term on the 
right-hand side of equation ( [TtI i is negligible. In this range one recovers the universal equation of the dynamical saddle-node 
bifurcation Q by taking X = xa~^^^ and T = ta}l^ , with the boundary condition Xit) « —\J—T at T tending to minus 
infinity. This property concerns the rectangular domain drawn on Figure ©c), where x is small, x ^ a^l"^, and t extends 
from —a^^l'^ \a t ^ a^^l"^, located before the abrupt increase. Therefore the time extension of this domain introduces the 
intermediate time scale, 

to^a-^l^, (19) 

long compared to unity (the adimensional time scale tcqk for the duration of a seismic rupture) and small compared to <b = a^^ , 
the average recurrence time between earthquakes. 

Let us prove that the short time is of order unity, by matching the solution ^(T) of the universal equation to the solution of 
equation ( l20b below, in the vicinity of the critical point tc{a) ~ a^^l'^ t^- Because X(r) behaves like ^-^^ before it diverges, 

it follows that the solution x(t) behaves as w a-^i^'^t -t "large" values of 5t ^ t ~ a~^^^tc before the critical time. Using 
5t in this development as time variable, x{5t) becomes of order one when 5t becomes of order one too. When this happens, the 
term at in equation ( fTTb is negligible, therefore the solution of this equation which can be matched with the solution near the 
bifurcation is the solution of the integrable equation 



(20) 
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with the asymptotic behavior for very large negative times x{5t) ^ ~Ji' ^^^^ equation shows that the time scale for the 
earthquake rupture is of order one, because it has no explicit dependence with respect to the small parameter a. This result 
is confirmed by the numerics: For a small we find that the rising time of x(t) close to tc{a) (defined as the half-width of the 
slope ^ solution of equation ( fTTb ) is toqk ^ 2.5, independent of a. From the observational point of view, the catastrophe 
takes place during this time toqk of order one, because the displacement is of order one then, compared to the displacement of 
order a}/'^ taking place during time a~^/'^ typical of the "universal" transition process. The two solutions match in the range 
1 -C (^<5t) <C a^^/'^. Supposing that the physical fast time scale for earthquakes is t"^^ ^ 10s, the intermediate time scale is 
^phys ^ ^_i/3^phys^ which is a fcw hours for a ~ 10"^. 



B. Precursor effects due to an external noise 

With a noise source added, the dynamical equation ([TJ becomes, 

da; 

+ at + eC{t). (21) 

at 

Actually the effective noise amplitude is not equal to e close to the saddle-node, but depends on the value of the parameter a. 
Indeed for \t\ < to , the cubic term in equation (l22l i is negligible, and the equation reduces to 

^=x^ + at + eCit). (22) 
dt 

which may be written on a form ^ ^ + T + i{a)C{t), by setting X = xa"^/^, T = ta^/^, and e(a) = ea"^/^. Therefore 
the effective noise is larger than e in the rectangular domain of figure (|2}c). 

Let us study the fluctuations of the solution x{t) of equation ( |22] |. For a small noise input, the solution may be expanded in 
power of e as above. At first order it gives 



dt 

whose solution is formally 



"^^^ [2xo(<) - ixl{t)]xi{t) + CW, (23) 



xi{t)= diC{i)exp[g{t)-g{i)], (24) 



where g(t) is the time integral of the second derivative of the potential 9{t) = Jf [2xo(u) — 3xq{u)]. The standard 

deviation (Tx-^ (t) has to be calculated numerically. We expect it to display the same behavior as for the cubic case in the whole 
domain where x{t) <C 1 , i.e. before the transition, and close to it, because the potential is cubic in this region. After the 
transition, we expect that the fluctuation decreases, because the solution without noise becomes quasi-steady. This is confirmed 
by the numerics: as for the cubic potential, the fluctuations strongly increase close to the critical time tc- With respect to time, 
the maximum of occurs at time tc{a) for small noise. Therefore the variance of the signal fluctuations cannot be used as 
a precursor for predicting the transition. Such a correlation between the standard deviation of the fluctuations and the sudden 
change of the solution has been reported recently IH where the GPS geodetic signal, which can be assimilated to our xo{t), is 
shown to be strongly correlated to the seismic signal (which we see as related to the fluctuations xi{t)). Note that when the 
noise increases, the growth of the fluctuations occurs earlier and earlier, their maximum progressively shifting before tc- This 
shift becomes visible only for an effective noise amplitude larger than unity, that is physically outside the range of noise values. 

Consider the case of small effective noise, where the correlation function and the spectrum of the fluctuations x{t) — xo{t) 
are well described by the correlation function and spectrum of xi{t), respectively. The calculation of these functions requires 
some care because the system is not in a statistically steady state. Therefore the spectral density of the fluctuations depends on 
time and the correlation function 

Tx,{t,T) = (a;i(t-T/2)xi(t + T/2)), (25) 
depends both on t and on r. A time dependent spectrum is formally defined by the (real) Wigner transform 



SxAtJ)= dTe-^"f^xi{t-T/2)xi{t + Tl2)), (26) 

J —OO 

that has to be modified for numerical appHcations, by introducing a sHpping window. 
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FIG. 3: (a) Width (in arbitrary units) of the correlation function of the fluctuation x{t) — xo{t) for a = 10 *^ {to = 10"^); (b) Spectral width 
(a. u.) 

The width of the correlation function and the spectral width A/ are reported in figures (|3) in a range of time of few to 
around the transition, together with the solution xo{t) drawn in solid red line for covering. Both widths show an interesting 
behavior which provides the same result. Consider first the left curve, obtained for the parameter value a = 10~^, typical for 
earthquake phenomena, as discussed above. The figure displays a strong increase of the correlation time of x{t) — xo{t) in 
the intermediate time range, reaching its maximum value at time t ^ 1.5 to (which was estimated as a couple of hours before 
the earthquake), then it displays a rapid decrease before the critical time tc{a) ^ 2340. In addition, we observe a slow growth 
of Tx as t increases from large negative values (not shown in the figure), increasing by a factor ten for — lOOto < t < 0, that 
corresponds to a time interval about one week. Such a remarkable behavior should be used as precursor The increase of the 
correlation length before the catastrophe can be understood when looking at the formal expression (l24l i. The second derivative 
of the potential vanishes at t = to, that leads to the flatness of g{t) in the whole domain < t < tc{a), as shown in figure (|4]-a). 
Figure (|4}b) shows that the "precursor time" (see caption) is given by the relation 



tprec ^ 2^0 



(27) 



which corresponds to about 4 hours. 
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FIG. 4: (a) Flat shape of g{t) before the critical time, for a = 10 ^, or to ~ 1000 ; (b) "Precursor time" tprec = tc(a 
Log scale ( ti/3 being the time where the width is i of its maximum value, before the catastrophe time tc{a)). 
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As for the spectral width, the result is just the opposite: it continuously decreases from large negative time, until the time 
t ^ l-5to, where it suddenly grows. The slow decrease of A/ corresponding to a slow shift of the spectrum towards low 
frequencies, is followed by a rapid spectral broadening at the end of the intermediate domain, before the transition time tc{a). 
The two stages of the width change are both important, because they occur before the transition. 

The growth of the fluctuations and their shift to lower frequencies can be understood as follows. As the transition approaches 
the potential V{x,t) becomes flatter and flatter, making weaker and weaker the restoring force toward the equilibrium. Therefore, 
at constant noise source, the amphtude of the fluctuations driven by this noise source will grow because the damping is less and 
less efficient. Moreover, the typical time scale for this damping will get larger and larger because of the decreasing stiffness of 
the potential, which will favour noise at lower and lower frequencies. 
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